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We establish a functional equation relating the Gibbs measure at a particular low temperature 
with the one at temperature 1. This equation enables us to calculate the limiting free energy of 
the Sherrington-Kirkpatrick spin glass model at this particular value of low temperature without 
making use of the replica method. We get additionally a sharper lower bound for the ground state 
energy e 

e > -0.7833 

close to the replica symmetry breaking and numerical simulations values. 
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INTRODUCTION AND MAIN RESULT 



During the last decade, mean field models of spin glasses have motivated increansingly many studies by physicists 
and mathematicians [3, @, H EH . The rigorous understanding of the infinite volume limit of thermodynamic 
quantities remained quite insufficient until the recent breakthrough obtained by Guerra and Toninelli [H on their 
existence and uniqueness. This major discovery followed by several important results 0, 0] providing a mathematical 
interpretation of the original formulae proposed by Parisi [8j on the basis of heuristic arguments. 

In this note, without making use of the replica approach, we calculate, for a particular value of the (low) temper- 
ature, the limiting free energy of the Sherrington-Kirkpatrick model and obtain a lower bound for the density of the 
ground state energy. Although the limiting free energy is given, for the whole low temperature region, by the rather 
complicated Parisi formula, we obtain, for a given temperature, a very simple expression. This allows an improvement 
of all known rigorous bounds for the ground state energy. 

We first recall some basic definitions. Suppose that a finite set of n sites is given. Let <n € {1,-1} be the spin 
variable on the site i and a a generic configuration in the configuration space £„ = { — 1,1}". The finite volume 
Hamiltonian of the model is given by the following real- valued function on S n 

Hn ^ = ^2 JijWji 

* l<i<j<n 

where the family of couplings Jij are independed centered Gaussian random variables of variance 1. Note however, 
that the function H n can equivalently be defined directly on the infinite configuration space Eoo = { — 1, 1} N and the 
infinite collection of random variables J = (Jy )i jS N just by the trivial modification in the definition of H n : 



v l<i<j<n 

This remark will be instrumental in the course of the proof since only changes the sequence of functions (i? n )n.eN and 
not the configuration and environment spaces. Henceforth, the Hamiltonian will be defined on these infinite spaces. 
For the inverse temperature (3 — y > 0, the disorder dependent partition function Z n (/3), is defined by 

Z n (0, J)=J2 eM-pH n (a, J)). 

(7 

Moreover, if Ej denotes the expectation with respect to the randomness J^, it is very simple to show that 
E J Z n (j3,J) = 2 n e^ n - 1 \ 

We denote by ^i n ^{a\J) 1 the corresponding Gibbs probability measure, conditionned on fixed randomness: 

finA^\J)=e- f3H ^ J) /Z n (/3,J), 

and, by S(fi nt p((x\J)), its entropy, defined by S(nn t p((r\J)) = - n n ,p(cr\J) log /J, n ,p((r\ J). 
For fixed randomness, the real functions 

f n {p) = —Ej log Z n (/3, J) 
n 

and 

f n (/3) = -logEjZ n (/3,J), 
n 

define the quenched average of the free energy per site and the annealed specific free energy respectively. 
The ground state energy density — e„(J) is defined by 

-e n (J) = - inf H n (a, J), 
n <re£„ 

For the low temperature region (/3 > 1), the J- almost sure existence of the infinite volume limits 

lim f n (fi) = fM, 
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and, 

- lim e„(J) = lim ^""(lO = -eo 

was first proved by Guerra and Toninelli [f| . More recently, Aizenman, Sims and Starr 0] gave a clear mathematical 
interpretation of the limit foo (/3) in terms of the variational formula proposed by Parisi. 
In the following section we prove the 

Theorem : Let (3* = 4 log 2 = 2.77258- • • Almost surely, the infinite volume limit foo(P*) is given by 

/ooGS.) = lim -E., log Z n ([3*,J) = ^ + ]- 

n-s-oo n 4 4 

A directly related result is the following corollary, which improves all the rigorous lower bounds for the ground 
state energy. 

Corollary : Almost surely, the ground state energy density of the Sherrington- Kirkpatrick spin glass model is bounded 
by 

e >= -0.7833 

PROOF OF THE MAIN RESULT 

Notice first that for all B> 0, the limit foo(P) exists and it is a convex function of ft jB[. Let /3 = fti = 1. From the 
high temperature results [l|, we have, almost surely, that 

/ oc (/3 1 )= lim -E J \ogZ n (/3 1 ,J) = f 00 (p 1 )= lim - log EjZ n (0x, J) = log 2 + ^ = log 2 + \. (1) 

The following figure [1] illustrates the definition of the inverse temperature /3*; the annealed free energy foo(P) — 
log 2 + is plotted as a function of j3 and the straight line is defined by 4^foo(Pi) = ft foo (Pi)- 

One can remark that, for j3„ — 2, 77258 • • • = 41og2, the annealed free energy foo(P*) is simply related to foo(Pi) 
via the following relation 

UP,) = f + log2 = + |ilog2) = |(log2 + \) = |/oc(/3x). (2) 

By making use of this remark we define the Gibbs probability measure fJ. n ,p,(o~\ J) by the functional equation 

cxp(p»ff»(q, J)) |* zj* (Pi, J) 



(3) 



induced by the mapping T : exp(/3 1 _ff„(cr, J)) (exp(/3 1 7J n (cr, J)))^*/ /Sl among Boltzmann factors. Since /i n>/ 3» is a 
probability on the configuration space, summing over a, we have indeed 

lim -Ejlog^^^J) = - aoo A), (4) 

where the limit a oa (^) gives the deviation of the free energy foc(P*) from its mean value: 

«oo :=ttco(^)= lim -^-—Ej log Z n (fii, J) — lim —Ej log Z n (P*, J) 

n->oo p! n n->oo ?i 

= ^foo(Pl) - /ooCS.) 

Pi 

= /oo(/?*) - /ooCS.)- 
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FIG. 1. The value /3* = 2, 77258 ■ • ■, is given by the intersection of the graph of the annealed free energy /oo(/3) with the straight 
line /?/oo(l). 



The reader can remark that the value /3i fixes the temperature scale; the temperature /3* is expressed in units where 
Pl = l, Therefore, ^ is a non dimensional quantity. 

We can now calculate the limit /oo(/3*)- The idea is the following. For r\ > 0, consider the set S of configurations a 
where the Boltzmann factor is close to its mean value: | exp(~f3H n (a, J)) — exp(^-(n — 1)| < rj. For /3 = 1, the set S is 
of full measure and moreover, the Gibbs measure behaves, at the thermodynamic limit, as linin^oo A log fi n ^ 1 (er| J) = 

?-/oo03i) = -log 2. 

Taking now the image under T, we have the following behaviour of the limit 

lim — log [i n ,i3, (&\J) = —7— -/oo(/3*) 

i.e. this image differs from the limit —fooifi*) by f 1 = log2. Now, from the previous equation and the easily verified 
relation 

^ "(«oo + f ))=/oo(/3*)-/oo(/3l), (5) 

we can deduce the behaviour of the limit by making use of the following argument. First, the difference of the 
limit — + from — log 2 + is given by (/3* — 1) log 2. Moreover, the difference between the limiting values 

- — + c*oo + — and - log 2 + aoo + x is also (/3* - 1) log 2. 

Second, from equation ([5]), since the limit — — + ctoo + — is equal to — /oo(/3*) + /co(/?i), it follows that — log 2 + 
ctoo + — = 0. We have indeed, 

10 ™ 

«oo + T =log2 = ^p, (6) 

and 

/ao(/3*)-/oo(/3i) = (/3*-l)log2 (7) 
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as illustrated in figure ([2]). 

One can now check that the obtained value of the limit 



o2 gl p2 

foo(M =P. log2 + ^ = !± + P± = 2.1718 ••• , 



(8) 



is lower than the spherical model bound (2.2058 • • •). 



lim„^ 00 -/oa(/3l) -log2 
< 



-/3,/4 + aoo = -/oo(/3.) + U>g2 



-«/4) 

-/oo09.) 
, .lim,,-,.^ i log (in^, (ct) 



FIG. 2. The straight oblique line represents the effect of the mapping T on measures. 



Remark 1: Our result can also be obtained using a large deviations approach, namely by comparing the Gibbs 
measures at the inverse temperatures 1 and /?* when the corresponding Boltzmann factors exp(— lH n (a, J)) and 

exp(— j3*H n (<j, J)) behave as their mean values exp( a ^i) and exp(^(n— 1)) respectively. In this case, the two limits 
lim„_>oo ilog/i ni| g,(o-|J) and lim ^log/i„ >/Sl ((j|J) differ by 

Remark 2: The point where these two limits are equal corresponds to the fixed point of the functional equation 
d3]) and arises when lim^oo i \og^i n ^ r (er| J) = - log 2 + . 

We can now obtain a lower bound for the ground state energy density — e„(J). Notice that 



fco (/?* 



lim — ) H n (a\J)fi n ,p,(a\J) + s(fip,) = -p 

n— >oo 71 * — • 
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1 

4' 



(9) 



where the limit 



s(M/3.)= lim —S{Hn,p, ( (7 I"0)> 

n— foo 77, 



(10) 



gives the (mean) entropy of the Gibbs measure. 
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We have indeed, by the positivity of the entropy s(//^, ) > 0, that 



eo = — lim 

P-HX3 



/oo(/3) 




(11) 



close to the value —0.7633 • • • obtained by numerical simulations based on the replica approach. 



CONCLUDING REMARKS 



In this note we obtained, under the assumption of minimal entropy, a rigorous lower bound for the ground state 
energy density which improves all the previous estimations. 

A last observation concerns the value of the temperature /?*: it is obtained from the relation ([2]) between the free 
energies /oo(/3*) and /oo(l); moreover, one can readily check that /3* is given by j3„ — (3%, where (3 C = 2y / Tog~2 is 
the critical temperature of the Random- Energy Model (REM) The REM is defined by 2™ energy levels Ei{i = 
1, • • • , n), a family of random, independent, identically distributed random variables; many results are qualitatively the 
same as those of the SK model. It would be interesting to clarify this relationship in order to obtain some information 
on the behaviour and properties of the Gibbs measure at low temperatures. Both f3 c and /3* are to be compared with 
the value at /5i = 1, i.e. the maximum value of (3 where the free energies of the two models coincide. What we learn 
by the comparison of the two models is that the Gibbs measure of the SK has seemingly a richer structure than for 
the REM. As a matter of fact, the entropy of the REM vanishes at f3 c while the entropy of the SK model is still 
strictly positive at this point. We expect moreover that the entropy of the SK model vanishes at /3* but this remains 
an open problem that is discussed in a forthcoming paper Q. 
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